Reconstruction of k-essence model 
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We explicitly construct the k-essence models which reproduce the arbitrary FRW cosmology, 
that is, the arbitrary time-development of the scale factor or the Hubble rate. The k-essence model 
includes scalar quintessence model, tachyon dark energy model, ghost condensation model as special 
cases. Explicit formulas of the reconstruction are given. 

First we consider the case that the action only contains the kinetic term. In this case, we find 
that the model reproducing the development of the universe in the exact ACDM model cannot be 
constructed although there is a model which reproduces the development infinitely closing to ACDM 
model. We find, however, the solution is not stable. 

Another is more general case including potential etc., where we find that there appear infinite 
number of arbitrary functions of the scalar fields, which are redundant to the time development of 
the scale factor. By adjusting one of the redundant functions, however, we can obtain the model 
where any solution we need could be stable. 

PACS numbers: 95.36.+X, 98.80.Cq, 04.50.Kd, ll.10.Kk, ll.25.-w 



I. INTRODUCTION 

Several observations tell that the expansion of the present universe is accelerating QJ-y] . In order to explain the 
acceleration, many kinds of models have been proposed. So-called k-essence model (j-Q is one of these models. The 
k-essence model is derived from k- inflation model 0, It is possible to reg ard the tachyon dark energy model [9rfl2l|. 
ghost condensation model (l7l H8|. and scalar field quintessence model |l3l4l"til| as variations of the k-essence model. 

In this paper we consider the reconstruction of k-essence model. We explicitly construct the k-essence models which 
reproduce the arbitrary FRW cosmology, that is, the arbitrary time-development of the scale factor or the Hubble 
rate. For general reconstruction, see [2l|, [22] . We consider two cases: One is the case that the action only contains 
the kinetic term and another is more general case including potential etc. In the former case, we find that the exact 
ACDM model cannot be constructed although there is a model infinitely closing to ACDM model. In the model, 
however, the solution corresponding to ACDM model is unfortunately not stable. In the latter case, we find that 
there appear infinite number of redundant functions of the scalar fields, which are not directly related with the time 
development of the scale factor. By adjusting one of them, however, we can obtain the model where the solution we 
need could be stable. At present, the roles of the other functions are not clear but we may obtain a model satisfying 
other constraints from cosmology by adjusting these functions. 

II. PURE KINETIC MODEL 



We may start with the following action: 

S = J d 4 x^g (A - K{q^)d^d^ , (1) 
where K and q are adequate functions. In the FRW background 

ds 2 = -dt 2 + a{tf ( dxl ) 2 > ( 2 ) 

«=1,2,3 

we may assume the scalar field cj> only depends on the time coordinate t. Since the redefinition of <f> can be absorbed 
into the redefinition of q(<fi), we further identify the scalar field <j) with the time coordinate t. Then we obtain the 
following FRW equation: 



^H 2 = K(-q(t)) + 2K'(-q(t))q(t), (3) 



and the equation given by the variation of <f>, we obtain 

= 2K"(-q(t))q(t)q'(t) - (6Hq(t) + q> ' (t))K> \-q(t)) . (4) 
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Here H is the Hubble rate defined by H = a(t)/a(t). Eq. Q can be integrated as 



1 ( Kr 



By differentiating @ with respect to t and using ([5]), we find 

G 



2 H = -6K a(t)- 6 ^q(t), (6) 
which gives 

(7) 

ft J\ Q 

Eq. ([7]) gives the explicit form of q(t) and can be solved with respect to t by using an adequate function / as 

t = /(?) • (8) 

Then ^ gives an explicit form of K'(—q) as 

^'Hz) = ^a(/(<7)r 3 - (9) 

Then for the arbitrary time-development of a or H given by a = a(t), Eqs. and @ give the forms of the functions 
K and q realizing the time-development. We should note, however, H cannot change its sign as we can see from (jSJ. 
Then the transition between non-phantom phase (H < ) and phantom phase ( H > ) does not occur. 

As an example, we consider the model corresponding to the Einstein gravity where the universe is filled with only 
one kind of perfect fluid whose equation of state (EoS) parameter w is constant. In the universe, the scale factor a 
and the Hubble rate H behaves as 



a(t)=a tm, H = ^^, (10) 



with a constant qq. Then Eq. (jJJ gives 



and Eq. (j9]) has the following form: 

K'(-q) = K Q q Q ^a- 3 q^- l i . (12) 

We can integrate ([12)1 as 

K{-q) =K X - ^^K q-^a Q 3 q^ , (13) 

where K\ is a constant of the integration. By using (jTTJ) , we find 

Aw 4 

Then we can check that Eq. (j3]) is satisfied if K\ = and we find the Einstein gravity where the universe filled with 
only one kind of perfect fluid can be reproduced by the model ^ if we choose q(t) and K(—q) by pd|) and (fT3"|) with 
K x = 0. 

We now extend the model (fTJ) to include the matters with constant EoS parameters Wi. Then since the energy 
density of the matters is given by poid~ 3 ( 1+Wi ^ with constants poi, the FRW equation ([3]) is modified as 

- 2 H 2 = K(-q(t)) + 2K'(-q(t))q(t) + £ p 0l a-^ 1+ ^ . (15) 



(i 
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Eq. (0| given by the variation of 4> is not changed and we obtain §5§ again. On the other hand, Eq. © is changed to 

6 



-H = -6K a(ty 3 y/q~(t) - J2 3 (! + Wi)p Qi a 



-3(H 



(16) 



which gives 



i(t) 6 (H'(t) + £ + m) P o t a-^ 1+w ^) 2 



(17) 



Eq. (jTTJ) gives the explicit form of q(t) as in ([7]) and can be solved with respect to t by using an adequate function / 
as in again. Then Eq.JS]) gives an explicit form of K'(—q) by (0). 
If we substitute the evolution of scale factor a(t), 



a(t) = Asinh 3 [at] , 



(18) 



which corresponds to the ACDM model in the Einstein gravity, to Eqs. (fl7|) and (fl2|) . then these equations give 
constant q and K. This result seems to be very trivial since K is nothing but the cosmological constant. Then we 
now consider the evolution of scale factor which differs a little bit from that of ACDM model in the Einstein gravity. 
For this purpose, we slightly modify (fTH| as 



a(t) = Aexp ln(sinh 3 [at]) + ef(t) 



(19) 



where A = (^f 1 ) 3 and a — |V3A are positive constants which reproduce the universe evolution in ACDM model, / 
is an arbitrary function of f, and e is a small constant. Then Eq. (JTTJ) gives 



«(<) 



^1 (^VL/ 2 - 3. 2 ^ sinh 2 H//" + ^ sinhVK/ 



//\2 



where we neglected the higher order terms of e. As an examples, we may consider the following function as /, 

/(f) = In (sinhl [af]) s.t. a(t) = A(smh% [at]) 1+t . 



By substituting ((3TJ) to ([20]). we find 



t ~ — sinh 1 



exp 



Then by further substituting (|2"21 into Eq. ([5]), we obtain 



if(-g) 



i^oA f 
/ dqq 

POm J 



exp 



K 
p0m e 



-2(l+e) 



eAe 1 



1 + e 



C — exp 



2Xo(l + e)y/q 



P0m e 



(20) 



(21) 



(22) 



(23) 



Here C is a constant of the integration. Then we find the existence of the functions q and K, which give a(t) very 
close to that of ACDM model. 

We now investigate the stability of the solution. So we now start with the two FRW equations and the field equation 
obtained from Eq. ([1]), 



—H 2 = K{-q4> 2 ) + 2if'(-#W , 
-L( 3 H 2 + 2H)= 1 -K(-q<j ) 2 ), 



dt 



(2a 3 q*<pK'(-q<j) 2 )) = 0. 



From Eqs. J24J) and (J25J), we obtain 



K'(-Q)Q 



H 



(24) 
(25) 
(26) 

(27) 
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Here Q = qcj) 2 . Inserting this equation to Eq. ([21))) , we find 



d 
dt 



(**>-»£)-<>, 



which can be integrated and gives 



(28) 



(29) 



Here i^o is a constant of the integration. We now write the form of a(t) as a(t) = aoe h ^ and consider the perturbation 
of the function h(t): h(t) = h (t) + Sh. Then from Eq. (I^U)) . we find the variation of Q is given by 



SQ 



*0 Jah, 



Ustih, 



By combining Eq. ([30]) and Eq. ([27|) . we obtain 



A'"(-Q )Qo = - 



h 



2h Q 8h 



h o 



K2 ^fa + ^t) 



(30) 



(31) 



Here Q is defined by substituting h(t) = h Q (t) into (gSJ): Q = age 6 ' l °( t )/i (i) 2 /^o K4 - Substituting Eqs. (J2ZJ), (I5D1) . 
and (|5T|) into the equation, ^hoSh — (^K'(-Qq) — 2K"(— Qq)Qq^8Q, which is given by the variation of Eq. we 
obtain a rather simple differential equation for 5h(t), 



6h- 3hoh ? + ' Ko 6h + 3ho8h = 0. 

hn 



(32) 



The above equation only contains ho and its derivatives. Then if we have ho, even if we do not know the explicit 
forms of q and K, we can find the stability of the solution. In case that solution ho is given by (|19p . by choosing 
a = A, 



Eq. ([32]) has the following form: 



h = ln(sinhs [at]) +ef(t) 



2a 2 

Sh - . , - Sh + 0(e) = . 



sinh [at] 



(33) 



(34) 



Since 



sinh 2 [at] 

model is not stable. 



> 0, the solution (|33|) is unstable. Then unfortunately the solution which behaves close to the ACDM 



By defining x = Sh and y — Sh, Eq. (|32j) can be rewritten as 



1 ( x 

dt\y 



3/i,Q 3ft,p/io+ h o 



(35) 



In order that the solution could be stable, the real part of all the eigenvalues of the 2x2 matrix in ([33)) should be 
negative: 



5iA ± < , A± 



3/ipfro + ho 
h 



± 



\ 



3h h Q + hp 
h 



12h Q 



(36) 



which requires ho > 0. This tells that the solution corresponding to non-phantom era, like quintessence era or matter 
dominance era, is not stable. 
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We now investigate the stability when we include the matters as in (fT5j) . Then the equations corresponding to 
- have the following form: 

^H 2 = K(-Q) + 2K'(-Q)Q + V PQl a~ 3{1+w > ) , (37) 

i 

^(3H 2 2H) = \k{-Q) - i J2 "W 3 ^) , (38) 

|(2a 3 Q5 A -(-Q)) = 0. (39) 

Then by considering the perturbation h —> ho + Sh, we find 

= Sh-bSh-cSh, (40) 
b ^ 3 ^(3fepfe + fee) - | ^(1 + ^ Wo ^-^e-^+^-feo 
^o + 1 Ei(l + ^ )po,a - 3(1+u, ' ) e-3d+-.)' l o 
K 2 ^ftg - 3E;(1 + ^)p 0i ao 3 C 1 +^) e -3(l+ Wi ) hD £ ( 6w ^ + (l + V) 



2 Mo + ! Ei(l + ^)poia - 3(1+Wi) e-3(i+-i)fto 

« 2 I Eg + + «; J )p 0i p 0j - a - 3(2+u,<+ "' j) e- 3 ( 2 +""+^) h ° 



(42) 



The above expression is rather complicated but in principle, we can check the stability by using the above equations. 
As long as we check numerically, the solution (1551 which behaves as almost ACDM solution in the Einstein gravity 
seems to be unstable. 

III. MORE GENERAL MODELS 

Since in the model ([1} , the solution which reproduces almost ACDM is unstable, we consider a more general model, 
whose action is given by 



S = J d A Xyf^ (J^ - K (0, X) + L mat ter) , A EE 8^8^ . 



(43) 



Then the FRW equations are given by 

j_#2 = 2X dK< ^ x X) - K (0, X) + p mattor , ~(2H + 3tf 2 ) = K (0, X) + p mattcr (i) . (44) 

As in (|15p . we consider the matters which have constant EoS parameters lUj. We now choose <j> = t, again, then we 
can rewrite the equations in (j44JI in the following form 

K(t,-1) = --L(2H + 3H 2 )-J2^a- 3{1+W >\ ggiM] = _L jj + * £ (1 + w . } pQ . a -s(.i^> t ) . (45 ) 

i X=— 1 j 

By using the appropriate function g(<f>), if we choose 

KfaX) = -^(2 5 "(0) + 3 5 '(0) 2 )-^^ Ol a- 3(1+ ^ ) e- 3 ( 1 +-)^) 

+ (X + 1) | ^ 5 "(0) + i £ (1 + Wi ) p 0i a^ 1+Wi) e-^+ w ^ \ + f] (A + 1)" ifW (0) , (46) 

we find the following solution for the FRW equations (|4"4"| . 

H = <?'(i) (a = a e ff(t) ) . (47) 



6 



Here K^(6) with n = 2, 3, • • • can be arbitrary functions. The case that iC^ n '(0)'s with n — 2, 3, • ■ ■ vanish was 
studied in fiol . [20| and the instability was investigated. 

We now investigate the stability of the equations without matter. From Eqs. ([33]), we can derive the following 
equation which does not contain the variable g", 



1 2 TT 2 00 

■i^Lx = -— + ^^(( n -l)X -n-l) X(X + l) n - 2 K^ (0) 



l + X 



H 2 H 2 



(48) 



where y = |= . Using Eq. (|4"8j) . we can rewrite dy/dN = (l/H)dy/dt (N is called e-foldings and the scale factor is 
given in terms of N as a oc e^) in the form which does not contain g: 



dy = „ Y 1 - y 2 ( </>_ 
dN l + X \ X 



2 oo 



+y - ip E v x ) (( n - l ) x - n - 1) + mx + 1)] (x + 1)"- 2 ^(")(0) 



(49) 



n=1 



We now consider the perturbation from a solution (j) — t by putting (f> = t + S(j) in (|49[) . Then we obtain 



dN 



-3 - 



g 12 dN \ 6g 



T^-^g") 



(50) 



If the quantity inside [ ] is negative, the fluctuation 5<f> becomes exponentially smaller with time and therefore the 
solution becomes stable. Note that the stability is determined only in terms of and does not depend on other 
(n ^ 2). Then if we choose properly, the solution corresponding to arbitrary development of the universe 
becomes stable. 

We now investigate the stability when we include the matter. Then the equation corresponding to Eq. (|4"5|) has the 
following form: 



1 -y 2 

3 —X = 

l+X 



H 2+ H 2 



( (n - 1)X - n - 1 ) X(X + l) n - 2 K^ (0) 



K X-l 

' IP 2(X + 1) 



Pmattcr 



2H 2 



Pmatter 



K X 

ipY+ 



i E p°* a o 



(51) 



Then we find 



dy 

dN 



3X 



y 



i-y 2 I 
i + x yx 

K 2 ( x-i 



2 00 

*y J2 [(0 + vX) ((« -i)x-n-i) + mx + 1)] (x + i)"- 2 k(")(0) 



2H 2 X \ X + l 

.2 



n=2 

yX) - <f> I Pmattcr 



2H 2 



yPmattcr 



(52) 



When we include the matter, the situation becomes a little bit complicated since not only H but the scale factor a 
appears in the equation. Then we need equation to describe the time development of o. By defining <5A for convenience 
as 



<5A = 3 2J(1 + Wi)p Qi a(t) 



-3(X+Wi) 



k 2 ( 5a 



(53) 



we obtain the following equations, 



d 

dN V 8X 



<p=t,H=g'W 



A B 
C D 



SX 



(54) 
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Here 



A = 



hl 8K (V A_ g » £ (1 + Wi )p 0i a(t)- 3il+Wi) , 



D 



C 



D 



24^( 2 ) 



8K^ 



6g> 



dN 



■ In ■ 



K 

W 2 



+ W4 )Po 4 aW^ 3(1+tUl) . 



(55) 



Generally, 2x2 matrix must have negative trace and positive determinant in order that two eigenvalues could be 
negative since the two eigenvalues are given by \{txM ± v^trM) 2 - 4(detM)}. So we just need to calculate the 
determinant and the trace of the matrix for investigating the stability of the fixed point <f> = t, H — g'(t). Since the 
expressions in (|55jl are very complicated, we consider the case that the solution is given by the behavior mimicing 
ACDM solution ([TBI in the Einstein gravity and the matter contents are given in the present universe. Then we find 



a(t) 



or z — ' 

i 



A sinh 3 [at] , 
2 

-a coth [at] , 
o 

2 f 

3 sinh 2 [at] 
4 



9 



x 1.86 x 1(T 



sinh 3 [at] 



(56) 
(57) 

(58) 
(59) 



where A = {pmo/PA)^ an d a = Ky/3pX/2. Note that in (f56|) . ([57|l . and ([58]) . we neglect the contribution from 
radiation, and in (|5"9")) . there only appears the contribution from radiation. Therefore the expressions in ([56]) - (fST?)) 
could be valid at least when t > 10 9 years. By using (|5"6"|) - (JSHJ), we find the following expressions of the determinant 
and trace of the matrix in Ea. (j54")) : 



tr 



det 



A B 

C D 

A B 

C D 



j«f.^( 2 )' _ 4 coshjat] 
a 3 3 sinh 3 [at] 



2 cosh [at] 2 coth [ a ^ (8f|if ( 2 ) + 



3 sinh 2 [at] 



(60) 



3 sinh [at] 



-27 



72 - 36 2 18 

cosh [at] 



-6 x 1.86 x 1CT 
1 



cosh [at] 

1 

cosh 4 [at] 
1 3 



sinh[a£| „„ , , \ k 2 t ^ 9 \, 
g— L +36tanh[at] -^K {2 >' 



- 12 x 1.86 x 10" 



1 



1 



sinh 3 [at] J a 



+3 



cosh 2 [at]sinh 3 [at] 2 cosh [at] sinh [at] 

1 



sinh [at] cosh" [at] 



+ 4 x 1.86 x 10" 



sinh 3 [at] 



(61) 



Note that 1 < cosh[at] < 2 and < sinh[at] < 1.7 in evolution of the universe, For example, if we consider the case 
that is constant, then the trace of the matrix is always negative and the determinant is positive when > 
since 72^^( 2 ) and 3/(sinh 2 [at] cosh 2 [at]) are dominant terms in Eq. lfoTj) . Therefore even if K^ 2 ' is constant, the 
fixed point solution mimicing ACDM solution in the Einstein gravity becomes stable as long as K^- 2 ' > 0. 



IV. TACHYON MODEL AS A GENERALIZED K-ESSENCE MODEL 



We now consider the tachyon model |9h12| as a variation of the generalized k-essence model. 
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The action of the tachyon model is given by 



Sbi = t 3 J d^-gViT)^/! + PJ(T)d^T. (62) 

Here l s is the string scale and r 3 is defined by using the string coupling constant g s as r 3 = l/(87r 3 <? s ^ 4 ). Since the 
scalar field T is dimensionless, we define a new field (j> by <j> = kT and redefine the functions V(T) and f(T) as 
V(T) = -v((f>)/T 3 and f(T) = u}{4>)k 2 /l 2 s . Then the action (JHSJ) can be rewritten as 



d*xy/-gv(<l>)yjl + u^d^t . (63) 



We now consider the FRW background ([2J). Since the redefinition of the scalar field <fi can be absorbed into the 
redefinition of ui(<f>), later we may identify with the cosmological time t. Then by expanding the square root in the 
action (|6"2"|) by the power of 1 + d^d^cj) ~ 0, we find, 

y/l+u(<t>)X = v / l-w(0)+w(0)(X + l) 

- ^^ + ^Mm ix+1) -w^ {x+1,2+0 ( (x+lf )> (64) 

where X = d^d^cf). Then by comparing Sbi in (|63|) with Eq.(|43|) and Eq. (|46|) . we can identify 



v(</>) = \ -\ (2g"(<t>) + Zg'{4>?) -J2 w iP°>"<> 



-3(l+v>i) „-3(l+u>i)ff(<« 



x S 



-^W a +EW (1+,U<) e- 3(1+tt " )flC *H , (65) 



= a Zzn^A „ > (66) 



(67) 



" 8 5 (W) + 3.g'(0) 2 ) + E i ^AHao 3(1+u " ) e- 3 ( 1 +«'O ff (*) 

If we consider the case that the action Sbi reproduces the ACDM behavior in the Einstein gravity, by using Eqs.(|56p- 
USSJ, we find 

K {2) {t) ~ j >0, (68) 

6* -4 sinh 2 [erf] +4 cosh 2 [at] sinh 2 [ai] + 1.86 x 10~ 4 sinh 3 [at] 



a 3 2 cosh[crf] 



KT ( 1 

9 sinh 3 \at] I -4 + 4 cosh 2 \at] + 1 .86 x 10" 4 \ \ 

\ sinh3 [at] J 

x ( -12 + 12 cosh 2 [at] + 1.86 x 1CT 4 \ ] <0. (69) 

\ sinh 3 [at] J 

Then we can find the stability of the fixed point corresponding to the ACDM model (TT8"| by substituting the equations 
m (ISgjl into Eq.dSHI) and (jl)Tj) . Since the expression is, however, very complicated, it is difficult to check the signs of the 
determinant and the trace analytically. Numerical calculation, however, tells that both of the determinant and the 
trace of the matrix are negative and therefore the fixed point solution corresponding to the ACDM model is unstable. 
So this model cannot reproduce the stable evolution of exact ACDM model. 
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V. SUMMARY 



In this paper, we consider the reconstruction of k-essence model which includes scalar quintessence model, tachyon 
dark energy model, ghost condensation model as special cases. Explicit formulas of the reconstruction were given. 
First we considered the case that the action only contains the kinetic term. In this case, we find that the model 
reproducing the development of the universe in the exact ACDM model cannot be constructed although there is a 
model which reproduce the development infinitely closing to that of ACDM model. We find, however, the solution is 
not stable. 

Another is more general case including potential etc., where we find that there appear infinite number of functions 
of the scalar fields, which are denoted by K^ n '(<f>) (n = 2, 3,4, • • • ) and redundant to the time development of the scale 
factor. Although K^ 2 '(cf>) is also redundant to the time development of the scale factor, we found that by adjusting 
K^ 2 \4>), we can obtain the model where any solution we need could be stable. At present, the roles of the other 
functions K^ n \<p) (n = 3, 4, ■ ■ ■ ) are not clear but we expect that we may obtain a model satisfying other constraints 
from cosmology by adjusting these functions. 
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